Spin-Hall effect in two-dimensional mesoscopic hole systems 
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The spin Hall effect in two dimensional hole systems is studied by using the four-terminal 
Landauer-Biittiker formula with the help of Green functions. We show that the heavy (light) hole 
spin Hall effect exists even when there are no correlations between the spin-up and -down heavy 
(light) holes and when the F-point degeneracy of the heavy hole and light hole bands is lifted due 
to the confinement or recovered by the strain. When only a heavy hole charge current without any 
spin polarization is injected from one lead, under right choice of lead voltages, one can get a pure 
heavy- (light-) hole spin current, combined with a possible impure light (heavy) hole spin current 
from two transverse leads. The spin Hall coefficients of both heavy and light holes depend on the 
Fermi energy, device size and the disorder strength. It is also shown that the spin Hall effect of 
two dimensional hole systems is much more robust than that of electron systems with the Rashba 
spin-orbit coupling and the spin Hall coefficients do not decrease with the system size but tend to 
some nonzero values when the disorder strength is smaller than some critical value. 

PACS numbers: 85.75.-d, 73.23.-b, 71.70.Ej, 72.25.-b 
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Recently spintronics has become one of the major fo- 
cuses of the condensed matter physics. 1 Generating and 
measuring spin current is still a primary task of this 
field. The spin Hall effect (SHE) is considered to be a 
hopeful way to produce the spin current. This effect, 
originally proposed as the extrinsic SHE as it requires 
spin-dependent impurity scattering, was first studied by 
Dyakonov and Perel' 2 in the early seventies and lately 
by Hirsch 3 and Zhang. 4 It further evolves into the in- 
trinsic SHE very recently as it appears in clean systems 
with the intrinsic spin-orbit (SO) coupling, such as two- 
dimensional electron systems (2DES) 5,6 ' 7 and bulk hole 
systems 8 without any impurity scattering. 

In 2DES, the intrinsic SO coupling comes from the 
Rashba 9 or Dresselhaus 10 terms which correlate the two 
spin states of electrons and lead to nonzero spin corre- 
lations (a^aic-cr)- Calculation based on Kubo formula 
found that the spin Hall coefficient (SHC) is a univer- 
sal value e/87r in macroscopic 2DES with the Rashba 
SO coupling. 5 The fate of the intrinsic SHE in the pres- 
ence of disorder is an important problem and has raised 
a lot of controversies. On one hand, Burkov et al. 
and Schliemann and Loss 12 pointed out that the SHE 
only survives at weak disorder scattering in macroscopic 
2DES. On the other hand, more analytical and numerical 
calculations 13 ' 14,15 ' 16 showed that the SHE vanishes even 
for any weak disorder. 

Situations are quite different in the mesoscopic sys- 
tems: using the four-terminal Landauer-Biittiker for- 
mula, Sheng et al. 17 and Nikolic et al. 18 showed that 
the SHC does not take a universal value in a finite meso- 
scopic 2DES, but depends on the magnitude of the SO 
coupling, the electron Fermi energy, and the disorder 
strength. Moreover, Sheng et al. showed that when the 
disorder is smaller than some critical values, the SHC 



does not decrease with system size but instead goes to 
some nonzero values. This robustness to the disorder in 
mesoscopic systems is suggested from the effect of the 
boundary effects (leads). 15 Therefore a pure spin current 
(no charge current associated with the spin current) can 
be induced by unpolarized charge current due to SHE in 
mesoscopic 2DES. 17 ' 18 ' 19 

The SHE in macroscopic hole systems has also been 
studied. Murakami et al. 8 predicted that the pure spin 
current can be obtained in p-type bulk semiconductors. 
They further suggested that this SHE comes from the 
Dirac magnetic mono-pole in momentum space because 
of the fourfold degeneracy at the T-point of the valence 
band. Then Murakami 20 pointed out that the vortex 
correction that was reported to kill the SHE in 2DES 14 
is identically zero in the Luttinger model so that the 
SHE can survive the impurity scattering. This is ver- 
ified by numerical calculation in bulk lattice Luttinger 
model very recently. 21 On the other hand, Schliemann 
and Loss 22 and Bernevig and Zhang 23 studied the SHE 
in 2D hole gases in which the HH and the LH bands 
are no longer degenerate at the T-point (no mono-pole in 
this case). However, they considered the case where there 
are still correlations between spin-up and -down HH's by 
introducing additionally the Rashba SO coupling in the 
hole systems. Therefore, whether the magnetic mono- 
pole is crucial for the SHE in hole systems is still un- 
known. 

We notice that the study of the SHE in mesoscopic 2D 
hole systems is still absent. There arc a lot of differences 
between 2DES and 2D hole systems: In 2DES, spin-up 
and -down electrons are always correlated via the Rashba 
or the Dresselhaus SO couplings. Nevertheless, in 2D 
hole systems the situation becomes much more compli- 
cated: For (001) quantum wells (QW's) with small well 
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width where only the lowest subband is important, un- 
less one adds an additional Rashba spin-orbit coupling, 
there are no direct or indirect spin correlations between 
spin-up and -down HH's (LH's), i.e., (a^ (Xk-o-) = with 
a = ±| for HH's and ±\ for LH's. The spin-up HH's 
(LH's) are only coupled with the spin-down LH's (HH's). 
Adding strain can change the relative positions of the HH 
and LH bands and the T-point degeneracy can be recov- 
ered. It is interesting to see whether there is still SHE 
in the absence of any correlations between spin-up and 
-down HH's (LH's) and it is also interesting to see the 
role of magnetic mono-pole to the SHE in the hole sys- 
tems. In this paper, we study these problems by using the 
four-terminal Landauer-Biittiker formula with the help of 
Green functions. 

We consider a 2D hole system which is consisted of 
a square conductor of width L attached with four ideal 
leads of width L/2 without any spin-orbit coupling (as 
illustrated in Fig. 1) in a (001) QW of width a. Due 
to the confinement of the QW, the momentum states 
along the growth direction (z-axis) are quantized. We 
take only the lowest subband for small well width. Then 
the Luttinger Hamiltonian in the momentum space reads 
with the matrix elements arranged in the order of J z = i , 
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where P ± Q = ( 7 i ± 7 2 )(P 2 + P 2 ) + E+, and R = 
-V3[-f2{P^ - Py) - 2ij 3 P x P y }. It is noted that the off- 
diagonal terms R and R^ mix the HH with the LH. In 
real space the Luttinger Hamiltonian can be written in 
the tight-binding version as: 
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where i and j denote the coordinates along the x- and 
y-axis; 71, 72 and 73 are the Luttinger coefficients, and 
mo/ (71 ±72) are the effective masses of the HH and the 
LH in the x-y plane with mo representing the free elec- 
tron mass; t = —h 2 /2moa,Q is the unit of energy with ao 
standing for the "lattice" constant, and (71 ±72)^ stands 
for the hopping energy; — (71 =p 2^ 2 ) 7^2- 1*| is the first 



subband energy in the z direction; ejj accounts for the 
spin-independent disorder, which is a random value in 
the range [-W/2, W/2}. Additionally 
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(3) 



is the strain Hamiltonian where ej 5 , is the strain-induced 

energy with e| 7^ e| . 25 Therefore, by adding strain, one 
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may either further increase the separation between the 
HH and the LH or reduce it to recover the T-point de- 
generacy. Moreover, it is seen from the Hamiltonian that 
there is not any direct or indirect spin flip between the 
spin-up and -down HH 's or between the spin-up and 
-down LH's in the Hamiltonian. 

u 



FIG. 1: Schematic view of a four-terminal junction of width L 
attached with four semi-infinite leads without any spin orbit 
coupling of width L/2. 

The spin dependent transmission coefficient from /1 
terminal with spin a to v terminal with spin a' is 
calculated using the Green function method 26 T°° = 

Tr[r°G™' R TZ'GZ'° A ], in which = i[£* - S^] rep- 
resents the self-energy function for the isolated ideal 
leads. 26 We choose the perfect ideal Ohmic contact be- 
tween the leads and the semiconductor. G££ R and 

Gl^ A , which can be obtained from G R(A ^ = {E - H c - 
, arc the retarded and advanced Green functions 
for the conductor, but with the effect from the leads in- 
cluded. Here Y, R ( A "> represents the sum of the retarded 
(advanced) self-energies of the four leads. 

We perform a numerical calculation when the leads 
are connected to isolated reservoirs at chemical poten- 
tials E + V,,, (//. = L, U, R and D), with the differ- 
ences between each other caused by the applied gate 
voltages. The particle current going through the lead 
\i with spin a can be obtained by the Landauer-Biittiker 
formula 27 i; = f E,^, [T#% - TgV v \. The spin 

current is defined as if = ^(1$ — I p. 2 ) for HH's and 

T 8,n = f ( 7 I - ^) for LH ' S in thc lcad M and the SHC 
for the pure spin current is defined as 



H(L) 



I^ L) /{Vl-V r ). 



(4) 



When thc charge current of the hole = e{l? K2 ' + 

Ip 2 *' 2 ^] = and thc spin current /jf^ 7^ for the lead 
\x, then is a pure spin current, otherwise an impure 

one. 
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FIG. 2: Particle currents and SHC v.s. the Fermi energy E 

when the pure HH spin current is obtained. L — 12. (a) 
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FIG. 3: Particle currents and SHC v.s. the Fermi energy E 
when the pure LH spin current is obtained. L — 12. (a) 

Solid curve: cr SH U ; Dotted curve: I£ (= I D 2 ); Dashed curve: 

_I I 3 _3 

I v 2 {= Ip). (b) Dotted curve: I^(= I D 2 )\ Dashed curve: 
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We drive a unit HH charge current without any spin 
polarization into the left lead L (if* = I and if =0). 
In order to get pure spin currents of HH in the upper and 
down leads (U and D), one needs to find a set of suit- 
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able combinations of V^, which lead to 
This can be obtained by choosing Vr = for conve- 
nience, Vu — Vd due to the symmetry and Vu/Vl 

)/(Ti+T 2 ) with T 1 =T I 



and T 2 = T 
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a pure HH spin current when 



T T l | 



2J-ij£ )/(Ti + T2) 7^ 0. This relation can be satisfied 
thanks to the phase shift provided by the last term of 
Hamiltonian (1) when holes hop from site (i,j) to site 
(i + l,j± 1). Similarly one can obtain the pure LH spin 
current. The main results of our calculation are summa- 
rized in Figs. 2-4. 

In Fig. 2 a pure HH spin current is generated in the 
U and D leads of a 12 x 12 square conductor made from 
the unstrained QW, where the T-point degeneracy of the 
HH and LH is lifted. In Fig. 2(a), all the particle currents 
in these two leads are plotted as functions of the Fermi 



energy. It is seen that Ifi 



to I? v 



D 



-II), leading 



= —I^jj. Therefore pure HH spin currents with 



opposite spin polarizations are obtained in the absence 
of any correlations between the spin-up and -down HH's 
and the T-point degeneracy. The SHC (jg H v of the upper 
lead is also plotted in the same panel. One finds that 
the SHC depends strongly on the Fermi energy. It is 
interesting to see that when the energy E is high enough, 
LH's produce an impure spin current as shown in Fig. 

2(b), where Ifi — I D 2 7^ —I u 2 = —Ifc and hence lf xV = 
IP D 7^ 0. One can also obtain a pure LH spin current 
combined with an impure HH spin current in both U 
and D leads when a unit HH charge current without any 
spin polarization injected into the left lead L, by using 
the condition Ifc v = ifc D = as shown in Fig. 3 for 

the same conductor. It is seen from Fig. 3(a) that Ijj — 

_i _i_ 1 

I D 2 = —Iu 2 = —I'd: which result in the pure LH spin 
currents with opposite spin polarizations in the U and D 
leads. The SHC of the LH again varies with the Fermi 



energy. Moreover, from Fig. 3(b) one has Ig 
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And one obtains 



ljj — —II) which lead to an impure HH spin current 
with Il v = ll D ± 0. 

Therefore, both HH and LH pure spin currents can be 
obtained through the suitable combination of the applied 
voltages. This provides us a unique way to manipulate 
the hole spin currents. Moreover, if one further allows 
only the HH charge current through the right lead R, 
then pure spin currents of HH and LH can be obtained 
at the same time in the U and D leads. 

In order to check the robustness of the SHC of the hole 
system, we plot in Fig. 4(a) the HH SHC of the U lead 
of a 12 x 12 conductor with and without strain versus 
the strength of disorder W over 5000 random disorder 
configurations. The strain is chosen to recover the V- 
point degeneracy of the HH and LH bands (and thus the 
magnetic mono-pole). For the strain free and strain ap- 
plied cases, E = 16|t| and 20|t| to ensure that the energy 
always sits around 1/4 of the HH bandwidth from the 
T-point. It is seen that for strain free case, the SHC is 
negative when W is small and then becomes positive. It 
reaches to the maximum near W — 30|t|, and then de- 
creases to zero when W is larger than 100|£| . When the 
T-point degeneracy is recovered by the strain, the SHC 
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FIG. 4: (a) SHC average over 5000 random disorder con- 
figurations versus the disorder strength. Solid curve: SHC 
of HH without strain when E = 16|t|; Dotted curve: with 
strain when E — 20|i|; Dashed curve: SHC of electron sys- 
tem. L — 12. (b) SHC averaged over 1000 random disorder 
configurations versus L without strain for W = 20|t|, 40|i|, 
80|t|, and 120|i|. (c) same with (b) but with strain. 



is killed by the same disorder strength, which shows that 
the degeneracy does not change the robustness of the 
SHE. It is noted that there is no correlations between 
the spin-up and -down HH's in both cases. We further 
compare the robustness of the SHC of 2D hole systems 
with 2DES. In Fig. 4(a) we also plot the SHC from the 
Rashba SO coupling of electrons in a conductor of in the 
same size as holes versus the disorder strength over 5000 
random disorder configurations. As the effective masses 
of electron and HH are different, we rescale the energy 
of electrons by 71 + 72 in the figure. The energy is again 
chosen to be at the 1/4 of the band width from the bot- 
tom. The Rashba SO coupling tg — Q.l\t\ (before the 
rescale) is taken to be the maximum value which can be 
obtained in the experiment. 28 It is seen that the SHC of 
electron system approaches to zero at W = 30\t\ (after 



rescale), much smaller than the value of HH's. There- 
fore the SHC of hole systems is much more robust than 
electron ones. This origins from the stronger intrinsic SO 
coupling of the hole systems. We further examine the size 
effect by calculating the SHC for the HH as a function 
of the system size L under different disorder strengths in 
the cases with and without strain. As seen in Figs. 4(b) 
and (c), similar to the case of electron systems with the 
Rashba SO coupling, 17 the SHC does not decrease with 
the size but goes to some nonzero value when the disor- 
der strength is smaller than some critical value. However, 
when the disorder is large enough, then the SHC goes to 
zero with the size. We have also calculated the SHC of 
LH's and come to the same conclusion. 

In summary, we have performed a mesoscopic investi- 
gation of the SHE for holes in four-terminal (001) QW's 
with a small well width so that only the lowest subband 
is populated, the correlations between the spin-up and 
-down HH's (or LH's) are totally absent and the T-point 
degeneracy between the HH and LH bands are lifted un- 
less a certain strain is applied. We find that the SHE 
still exists. Moreover a pure HH (LH) spin current can 
be generated combined with a possible impure LH (HH) 
spin current, when a HH charge current without any spin 
polarization is injected from one lead. The SHC's for 
both HH and LH depend on the Fermi energy, device 
size and the disorder strength. We also find that the 
SHE of holes does not need the mono-pole from the T- 
point degeneracy. We show that the SHE of 2D hole 
systems are much more robust than that in 2DES. This 
is consistent with the latest study in bulk systems. 21 The 
SHC does not decrease with the system size but tends to 
some nonzero value when the disorder strength is smaller 
than some critical value, similar to the electron case but 
where there are direct correlations between the spin-up 
and -down electrons. 
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